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The long-standing problem of the asymmetry between matter and antimatter in the Universe is, 
in this paper, analysed in the context of the modified theories of gravity. In particular we study 
two models of f(R ) theories of gravitation that, with the opportune choice of the free parameters, 
introduce a little perturbation to the scale factor of the Universe in the radiation dominated (RD) 
phase predicted by general relativity (GR), i.e., a(t) ~ t 1 ^ 2 . This little perturbation generates a 
Ricci scalar different by zero, i.e., R ^ 0 that reproduces the correct magnitude for the asymmetry 
factor r) computed in the frame of the theories of the gravitational baryogenesis and gravitational 
leptogenesis. The opportune choice of the free parameters is discussed in order to obtain results 
coherent with experimental data. Furthermore, the form of the potential V, for the scalar-tensor 
theory conformally equivalent to the f(R) theory which reproduces the right asymmetry factor, is 
here obtained. 
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I. INTRODUCTION 


Observational data suggest that our Universe is com¬ 
posed for the most part of matter, while the antimat¬ 
ter is only presents in trace amounts [lj. Recent stud¬ 
ies propose that the origin of this asymmetry between 
matter and antimatter lies in the beginning phase of the 
Universe, before of the Big Bang Nucleosynthesis (BBN) 
Different theories introduce different interactions 
Beyond the Standard Model (BSM) in order to explain 
the origin of this asymmetry in the primordial Universe 
GH3- In this paper we will show how a f(R) theory and 
a particular non minimal coupling between Ricci scalar 
and matter are able to explain the cause of this asymme¬ 
try. In particular we have shown how a small correction 
to the standard Hilbert-Einstein action allows the repro¬ 
duction of a small variation of the scale factor of the 
Universe in order to reproduce the expected asymmetry 
factor. We will describe the ln(R ) gravities never intro¬ 
duced before in this context which represent a suitable 
alternative theory in order to describe the early Universe 
phenomenology. This work is organised as follows. In 
section II we will briefly introduce some important pa¬ 
rameters related to baryogenesis and to leptogenesis. In 
section III we will briefly resume the main topic of metric 
f(R) theories of gravity and their implications for Uni¬ 
verse dynamics. In section IV we will show how modified 
gravities can reproduce the correct baryon asymmetry 
factor by means of two functional form of f{R) proposed 
for the first time in this context. In section V we will 
introduce new results about leptogenesis for an already 
studied form of f(R) and for a new one, i.e. (IIV.31I) . 
originally proposed (in this context) for the first time in 
this paper. In section VI we will obtain the functional 
forms of the potential of a primordial scalar field which, 
in a scalar-tensor theory of gravity, can realize the lepton 
asymmetry of the same magnitude of that one generated 
by the f(R) analysed here. The technique, adopted in or¬ 
der to obtain this potential form, is based on the confor¬ 


mal equivalence between f{R) theories and scalar-tensor 
ones. Free parameters of this potential are fixed by fixing 
free parameters of the f(R) theories in order to obtain 
the expected asymmetry factor. In Section VII we will 
summarize and comment on our results. 


II. MATTER-ANTIMATTER ASYMMETRY 


In the longstanding list of attempts to explain matter- 
antimatter asymmetry in the Universe, several parame¬ 
ters have been introduced. An important parameter used 
to quantify the amount of baryon matter that exceeds 
antibaryon matter is the asymmetry factor r]B , i.e. the 
Baryon Asymmetry Factor (BAF): 


r]B 


n B - ng 


(HI) 


where rising) is the number of baryons (antibaryons)per 
volume unity and s the entropy density for the Uni¬ 
verse. Some works about Cosmic Microwave Background 
(CMB) anisotropies and BBN show that this factor is 
?7 ~ 10~ 10 fl]. Analogously, it is possible to introduce 
the ?y factor for leptons, i.e. the Lepton Asymmetry Fac¬ 
tor (LAF): 


Vl = 



(II.2) 


where n^ng) is the number of leptons (antileptons) per 
volume unity and s the entropy density for the Universe 
[lj]. For rjL there are not experimental constraints, but 
only deductions that estimate it with the same magni¬ 
tude of r]s [IJ- Another useful quantity is the baryons to 
photons ratio: 


n s /n 7 « 6* 1CT 11 , (II.3) 


or the ratio between quarks and antiquarks in the pri¬ 
mordial Universe (t < 10“ 6 s): 


n ^ n ^ 3*10- 8 , 

n q 


(II.4) 
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where n q , n q and n 7 are respectively the density of 
quarks, antiquarks and photons in the primordial Uni¬ 
verse. 


III. f(R) THEORIES OF GRAVITY 


Since the discovery of the current accelerated phase 
of the Universe Q, and the hypothesis of the early time 
inflation [j], many alternative models of classical or quan¬ 
tum gravity have been proposed. f(R) theories are one 
of the most significant attempts to explain the current 
expansion of the Universe, the Dark Matter behaviour, 
and inflation 0 ®. In this paper we will show how f(R) 
theories can be responsible for the asymmetry between 
matter and antimatter. f(R) theories are obtained by re¬ 
placing the Ricci scalar R in terms of a generic function 
f(R ), modifying correspondingly the Einstein-Hilbert ac¬ 
tion S as 

s = ^ Jd 4 x^gf(R) + S M , (III.5) 


where k = jp-, Mp is the Planck mass, g the determi¬ 
nant of the metric tensor involved and Sm the action for 
matter terms. 

In order to obtain dynamical equations, it is possible 
to vary with respect to the metric g the Eq. (IIII. 51) . 
obtaining the fourth order field equations 1,0 

f'(R)Rp V - \f(R)g - [V M V„ - g^D] f'(R) = kT^, 

(III.6) 

where 


= 


—2 SSm 

yf=g 5g» v 


(III-7) 


or equivalently, splitting Ure matter counterpart from 
curvature contribute [10, 11], i.e. G a p = R a p — ^ Rg a p = 

Tjpur-v) w j iere as shown in mm we define 

T { :r ] = j^{g«0im-Rnm/2+ 


+ f(Ry^ (g a ^gp v - g aP gp V )} (III.8) 

as the curvature energy-momentum tensor, where de¬ 
notes the covariant derivative with respect to the indices 
g and v. In this work we consider the Friedmann metric 
for the Cosmos: 

ds 2 = — dt 2 + a 2 (t) f dr 2 + r 2 dfl 2 J , (III.9) 


where a{t) denotes the scale factor of the Universe, k is 
the curvature of the space , and dfl 2 = d0 2 + sin 2 9d(t> 2 . 
Thus it is possible to obtain the modified Friedmann 
equations, i.e., II El 


rr2 & 1 

H + — =o 

a z S 




(III.10) 


and 


-2H-3H 2 ~\ = P curv + (III.ll) 

a 2 }'{R ) 

where the dot, i.e. ', denotes the derivative with respect 
to the cosmic time, and P m and p m are respectively the 
pressure and the density of all fluids which fill the Uni¬ 
verse. Besides the curvature density is defined as 

Pcurv = y lf(R) - Rf'(R)} ~ 3HRf"(R)^j , 

(III.12) 

and the barotropic Pressure is denoted by 

Pcurv — ^curvPcurv , (III.13) 

where the effective curvature barotropic factor is given 

by 


U) r 


= -1 + 


Rf"(R) + R Rf"'{R) - Hf"{R) 


[f(R) - Rf'(R)} /2 - 3 HRf"{R) 

(III. 14) 

In the following we will use the positive signature 
(—,+,+,+) and the Ricci scalar R will be written in 
function of the Hubble parameter as ns® 


R = 6 


H + 2H 2 


(III.15) 


According to PLANCK result EH, in the following we 
will impose the spatial curvature equal to 0, i.e. k = 0. 
If we denote the right side of Eq. (IIII.10I) and (IIII.11I) 
respectively as Peff — P (R) Pcurv and P e jj — ~f\R) 

P e ff , it is possible to introduce the EoS effective param¬ 
eter as: 


W e// 


Peff 1 2 H 

Peff 3 


(III.16) 


IV. BARYOGENESIS 

In 1967 Sakharov inferred three necessary conditions 
to obtain a net baryon asymmetry 0,®. 

These three conditions are: 

• existence of reactions violating baryon number; 

• violation of the C and CP symmetry; 

• the Universe needs to be out of the thermal equi¬ 
librium for a finite period of time. 

Later studies have shown how it is possible to explain the 
asymmetry relaxing some of these three conditions. For 
example in 1987 Cohen and Kaplan 0 proposed a model 
of spontaneous baryogenesis in which the CP violation 
and out of equilibrium phase were relaxed and a dynamic 
violation of CPT symmetry is introduced. Accordingly, 
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the expanding Universe breaks CPT symmetry, that is 
restored considering a static Universe. In this model an 
interaction between a scalar field <f> (called ilion and it 
is an axion-like particle) and the baryon number current 
Jg is introduced, i.e., 


jdM), (IV.17) 

where / is an energy scale and it is usually bigger than 
10 13 GeV. This term dynamically violates CPT in an 
expanding Universe. Considering Noether Theorem and 
considering <j> field uniform in space component we can 
rewrite (|IV.17I) as 



~ n B ). 



(IV.18) 


As explained in [3], this interaction generates a dif¬ 
ferent population of barions with respect to antibarions 
before the field </> reaches the minimum point of its po¬ 
tential. In this way the g factor is equal to 


V ~ 


g*fT ’ 


(IV.19) 


where T is the temperature at which the asymmetry and 
g*, i.e. the total number of freedom grades of all the 
fields present in the early Universe, are computed. No¬ 
tice that T is also the temperature to which the reac¬ 
tions, which violate the baryon number, decouple from 
the background. 


A. Gravitational baryogenesis for f(R) = R + aR 2 


Gravitational baryogenesis was first proposed by 
J.Paul Steinhardt et al.in 2004 jT8j. Inspired from spon¬ 
taneous baryogenesis jBj they proposed this kind of inter¬ 
action: 

(IV.20) 

where R is the Ricci scalar, Jg is the baryon current, 
M* the cut-off scale of the effective theory. This interac¬ 
tion emerges in the phenomenolo gy o f some theories of 
quantum gravity or supergravity I18 h 23I |. Relaxing the 
equilibrium condition it is possible to reproduce the cor¬ 
rect asymmetry because interaction (IIV.20I) dynamically 
violates CPT. In this case rj becomes: 


R 

V ~ MITd ’ 


(IV.21) 


where To is the temperature at which the decoupling 
of reactions that violate baryon number from primordial 


plasma occur. Tracing general relativity equations p4| . 
i.e. Rfj, u — bg^R = we can express Ricci scalar 

V1 P 

as: 


R = 8n 


(1 — 3 w)p 

Ml 


(IV.22) 


The baryogenesis happened in a radiation dominated 
phase, when the EoS parameter is 1/3, i.e. w = 1/3, 
so the Ricci and its derivative are equal to zero. If we 
hypothesize perturbative effects in QFT 0 , or a mod¬ 
ified gravity we can obtain an adiabatic index slightly 
different from 1/3. 

In the following, we will show how an f(R) theory can 
reproduce the correct value of the Ricci scalar and its 
first derivative, in order to explain the expected asym¬ 
metry factor. Indeed, we will search for an f(R) theory 
of gravity that can reproduce a scale factor slightly dif¬ 
ferent from the expected one for the radiation dominated 
phase, i.e., 

a(t) = a 0 t 1/2 + A (t), (IV.23) 


with A (t) little as much as to not change the known ther¬ 
mal history of the Universe. 

In the following, inspired by the work [25], [0 , which 
first suggested the idea of applying f(R) in order to ex¬ 
plain baryogenesis, we will study a f (R) form which can, 
differently from the one proposed in 0 , evade the solar 
system test 0, 01 , and that is more suitable to describe 
the early Universe dynamics, i.e., 

f(R) = R + aR 2 , (IV.24) 

where a is a constant having the dimension of GeV~ 2 . In 
particular, immediately hereafter we will show how this 
f(R) can generate a little shift of the scale factor, from 
the standard one, i.e. a = dot 1 / 2 in Radiation Dominated 
(RD) phase, where ao is a dimensional constant. 

We are looking for a solution of the scale factor given by 
Eq. (IIV.23I) , and in order to find this form of solution, we 
substitute the (IIV.23j) expression in the Friedmann Equa¬ 
tions (IIII.10I IIII. Ill) . We solve Friedmann Equations in 
the linear approximation for A (t). This linearization is 
performed because we suppose that A(f) is a little per¬ 
turbation with respect to the scale factor a(t). Indeed 
we will check the time lapse in which this linearization is 
true by the means of the introduction of the e parameter, 
i.e. 

e = - = e °(xr) 7 ’ ( IV - 25 ) 

a Mp 

with e 0 = (f.^r/ 2 X ° A ^ —. In the previous equa¬ 
tion we have used the relationship between cosmic time 
and temperature in the radiation dominated phase, i.e., 
t = { 32 + g , ) 1 ^ 2 ^ (with g * ft: 106). The linear approx¬ 
imation is valid if e <C 1. Besides we solve Friedmann 
Equations imposing w = 1/3 + 5 and the following ansatz 
for A (t): 


A(f) = A 0 


(IV.26) 
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FIG. 1: Baryon Asymmetry Factor (BAF) in function of 7 . 



FIG. 2: BAF in fucntion of £0 for f(R) = R + aR 2 . 


with Aq dimensional constant. 


In this model the Ricci scalar computed from the mod¬ 
ified scale factor a, always in the linear approximation, 
is equal to 

Q \ 

R= 2 ^ (4/32_1) ’ (IV ' 27) 

and its time derivative is: 

I?=^(4/3 1 2 -l)(/3-!)A (IV.28) 

In this way we can substitute this expression in the Eq. 
(IIV.20I) and using the relationship between time and tem¬ 
perature in the radiation-dominated phase, the parame¬ 
ter e defined in (1IV.25I) and imposing M* «Mpwe obtain 
the asymmetry factor as: 


V 


3 / 167r 3 g* 

4 £ ° \ 45 


3/2 

7(2 - 7)(7 + 4) 



(IV.29) 

where 7 = 1 — 2/3 and Tjj is the temperature at which 
the reactions which violate baryon number decouple from 
the background. Please note that eo is a dimensionless 
constant because of the different dimensionality of the 
constant ao and Xq. Given the arbitrariness of the ratio 
Ao/ao we can choose it as small as necessary in order to 
obtain eo equal to 1 as it will be done in the following. 
In Figure [I] and [2] we report respectively 77 in function 
of 7 and eo, for a determinate value of decoupling tem¬ 
perature, i.e. Tjj = 10 16 GeV. The model (IIV.24I) is 
able to reproduce a correct asymmetry factor for a wide 
range of 7 , eo- Furthermore the model consider the nec¬ 
essary condition e << 1. If we suppose £0 equal to 1 and 
Td — 10 16 GeV we get the results presented below: 


r 

1. xKT 10 - 
8. x 10~ n - 
6.XKT 11 - 
4. x 10~ n - 

2. XKT 11 - 


FIG. 3: BAF r) in function of the Temperature for S = 0.001 
(green line) for the model f(R) = R + aR 2 . Expected value 
for 7 is in yellow. 



1. r) ~ 1CT 10 , 

2 . e < 1 , 

3. 5 < 1/3. 


e = 10" 37 


0.5 per 7 = 0.1 
0.06 per 7 = 0.4 


(IV.30) 


In Figure [3] we show rj in function of Tp>. Our model 
is consistent with observational cosmic history because it 
verifies: 
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B. Gravitational baryogenesis in ln(i?) theories 

In this section we describe a f{R) model, originally 
proposed in this paper as a suitable attempt to reproduce 
the expected BAF. In particular this modified gravity can 
produce, in a very simple way, a slight variation of the 
scale factor with respect to the scalar factor of the radi¬ 
ation dominated phase predicted by General Relativity 
(GR). We consider again the presence of the interaction 
(1IV.20I) that allows the splitting of the energetic level be¬ 
tween matter and antimatter. The f(R) model in exam 
is Q: 


f{R) = R + 1 R~ n {ln^) m , (IV.31) 

where n is a real number, i.e. n > — 1 , m is arbitrary, 7 
and g are some dimensional constants. It is easy to show 
that the solution of Friedmann Equation for this f(R) is 
0 : 


Q + l)(2n + l) 

a^t "+ 2 


(IV.32) 


with the effective EoS parameter (defined in (1111.161) 1 
equal to: 


v 



FIG. 4: BAF in function of n. For n of the order of 10 5 
(a ~ t 0 50001 ; we reproduce expected BAF. 


We obtain results with the expected BAF, adopting n 
smaller than 0.01, which reproduces a scale factor like 
£ 0.5001 or £0.50001 ^ j n p er f eC £ agreement with the thermal 
history of the Universe. In Figure [I] we show the be¬ 
haviour of 77 in function of n, for some fixed values of 
other parameters. 


6 ?r 2 + 7n — 1 

Weff = ~3(n+l)(2n + l)' 


(IV.33) 


V. GRAVITATIONAL LEPTOGENESIS 


It is easy to show that for n = 0 scale the factor has 
the same behaviour of GR scale factor in the RD phase. 
Straightforwardly the power-law of the solution depends 
only on n and not on m and other two dimensional con¬ 
stants. 

Substituting (IIV.32I) in (IIII.15I) , and deriving the result 
with respect to the cosmic time, we can get R , i.e., 


• 12n(5 + 19?r + 22n 2 + 8 n 3 ) 

~ (2 + n) 2 t 3 


(IV.34) 


which is equal to 0 for n = 0, as expected. Rewriting 
(IIV.34I) in terms of the temperature, we can rewrite all 
as: 


12n(5 + 19n + 22?i 2 + 8 n 3 )T 6 
(2 + n) 2 0.027<?*“ 3/2 Mj, 


(IV.35) 


In this way the asymmetry factor (IIV.21I) becomes: 


12n(5 + 19n + 22n 2 + 8 n 3 )T 5 
(2 + n) 2 0.027g* 3/2 MpM% 


t d 


(IV.36) 


For example, if we choose n = 0.01 that reproduces a 
scale factor proportional to t 051 and impose decoupling 
temperature equal to 10 16 GeV, we get the following 77 
factor 


77 = 


6 • 10 26 


GeV 2 . 


(IV.37) 


If we suppose that the Universe has a baryon asym¬ 
metry, it is straightforward to speculate about a lepton 
asymmetry. Common sense leads us to hypothesize that 
this asymmetry is of the same order of magnitude as the 
baryon one. In particular the charge neutrality of the 
Universe is a good evidence to support this assumption, 
even if one of the main problems is that we do not have 
any direct measurement about the lepton number given 
by three neutrinos. Predictions of BBN assume that neu¬ 
trino lepton number is very small [T[. 

In the last few decades different models of leptogenesis 
have been developed l!, :2()l - (35t . In the following we 

will consider a theory that explains lepton asymmetry as 
a consequence of scalar curvature, similarly as we have 
proceeded for the gravitational baryogenesis [Hi, [H, [35}. 
In the model [31!, [H, [HJ the origin of lepton asymmetry 
is realized by a new lagrangian in which a coupling be¬ 
tween Majorana Neutrinos and Ricci scalar. This term 
violates CP is present. Hereafter we will explain the basic 
idea behind this interaction. We start realizing that in 
the Standard Model only neutrinos with defined chirality 
exist, i.e., left-handed neutrino (LHN) and right-handed 
antineutrino (RHA) . Some extensions of the Standard 
Model introduce heavy right-handed neutrino (HRN) and 
heavy left-handed antineutrino (HLA). These neutrinos, 
whose masses are described by the see-saw mechanism 
[13, [37], are main actors in the model of leptogenesis stud¬ 
ied hereafter. The above introduced interaction, which 
violates CP, and that it is even under C, and odd under 
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P, is described by [M: 

Ay = V-gOR-ipi isi>, (V.38) 

where 0 is a constant that introduces the effective range 
of validity of the model, i.e. Mp 1 , ip is a fermionic field 
(neutrino in our case), i is the imaginary unit, 75 the 
chirality operator . 1 

The operator (IV.381) conserves CPT only in a static Uni¬ 
verse and dynamically violates it in an expanding Uni¬ 
verse, i.e. R / 0. 

Introducing (IV.38I) the energy of a fermion is determined 
by its chirality through the interaction with the gravi¬ 
tational background. Euler-Lagrange equations give this 
result for dynamics of fermions coupled non minimally to 
the background 2 


i'y^dn'ip — Mip — iOR^ip = 0. (V.39) 

Energetic dispersion relation is given by (see |3l], [35} 
for clarifications): 

E 2 ip = (p 2 + M 2 + 0 2 R 2 )ip - 9{'y B 'fd ll R)ip. (V.40) 

Now we express ip as a superposition of a left-handed 
spinor ip- and of a right-handed one ip + , i.e. 

ip = ip-+ip + , (V.41) 

where 

f 1 —I— 'y 3 * ) 1 — q 3 

iP + = ^-M-R-ip, ip _ = — Y~ip. (V.42) 

We use this property for 75 : 

(£) = (-£)• <™> 

It is easy to show that, through chiral properties (IV.43I) . 
the effective energetic levels of spinors are split in func¬ 
tion of their chirality [31, 35), i.e., 

E 2 ip± = (p 2 + M 2 + d 2 R 2 ± 6R)ip±. (V.44) 

From (IV.44I) we get this expression: 


E± = Vp 2 + M 2 + 9 2 R 2 + 


OR 

2a/p 2 + M 2 + 0 2 R 2 ' 

(V.45) 


1 The 75 is the product of four Dirac matrices: 75 = 270717273 - 
Trough Wick rotation we can recast it in the Euclidean space 
as 75 = 70717273 - Dirac Matrices allow to build fermionic 
fields observables and they satisfy anticommutation relationship 
{7m 5 7 ^} — 2 Ig^iy. 

2 The equation is obtained by varying the action composed by the 

Dirac term plus the term which violates CP, i.e: C = 7 ^ — 

+ 0 R'ipi'ys'il). 


If we consider Major ana neutrinos in the chiral form 
(Np,Np) T we can denote the first component of the 
bispinor Nr as the HRN while the second one as the 
HLA, i.e. N^ = N l . 

If we denote Nr with ip + and Nl with ip- the energetic 
levels of this particle are given by (IV.441) . In particular 
we point out that effective masses of two neutrinos are 
different by means of the gravitational interaction. 

Hence from (IV. 441) and for OR C M we get: 


M±=M + 


o 2 r 2 or 

2 M ± 2 M ’ 


(V.46) 


that shows how left-handed fermion component has 
an effective minor mass respect to the right-handed one. 
Furthermore the rate of decay for a massive heavy neu¬ 
trino is S3, HI: 


r± = - l -h 2 M ± , (V.47) 

where ± denotes particle chirality and h is Youkawa cou¬ 
pling. 

Lepton asymmetry factor can be witten as 



r+ - E_ M+ -M- R 

r+ + r_ “ m ~ 


(V.48) 


in the limit of OR « M. In the following we will con¬ 
sider 0 « Mp 1 . 

Please note that in (IV.481) the mass of the heavy neutrino 
is present. This mass, nowadays, has not a certain value. 
Different theoretical models proposed a wide range for 
this mass that goes from 10 GeV to 10 16 GeV. Here¬ 
after we will estimate the mass of the heavy neutrino in 
function of the decoupling temperature, i.e., TA, for the 
lepton number violation reaction. 


A. f{R) = R + aR 2 applied to leptogenesis 


In the following, we will compute the correct lepton 
asymmetry factor (LAF) in f(R) = R + aR 2 gravity and 
we will give a probant prediction of the mass of the right- 
handed neutrino. If we substitute (IIII.15I) in (IV.481) and 
we express everything in function of temperature we get 
the following expression for LAF: 


0=^0 


( 167r 3 g» 
V 45 




M 

(V.49) 

In Figure[5]we reproduce i] in function of the mass M for 
£o = 10 _9 ,7 = —1, Td = 10 15 GeV. The correct LAF 


rj eX p ~ 10 ~ 10 is obtained for M ~ 10 12 GeV. 

If we increment the Decoupling temperature (till to 
10 16 GeV) we see that the mass of the heavy neutrino 
necessary to produce the correct LAF without violate 
BBN constraints [35| is of the order of 10 16 GeV. 
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H 



M(GeV) 


FIG. 5: LAF r/ in function of HRN for f(R) = R + aR 2 . 


Now through see-saw mechanism |36,, [37] we estimate 
light neutrino mass from heavy neutrino mass. The rela¬ 
tionship between light and heavy neutrino from the see¬ 
saw is: 


m v M ~ m 2 D , (V.50) 

where rajj is the Dirac Mass, i.e., the mass of one of 
standard model fermions. Considering our best fit for 
neutrino equal to M = 10 12 GeV and choosing for tod 
quark bottom mass (4 GeV), we get our prediction for 
light neutrino mass: m v ~ 10 _1 eV a value in perfect 
agreement with experimental constraints from neutrino 
oscillation in the atmosphere [37}- Analogously, if we 
apply the same mechanism to our second fit for heavy 
neutrino (M = 10 15 GeV), and choosing as tod quark 
top mass ( 10 2 GeV), we get another consistent value for 
our prediction of light neutrino mass (to,, ~ 10~ 2 eV). 

At the end we should emphasize a process in which 
it is possible to obtain baryogenesis via leptogenesis, or 
rather a baryon asymmetry from a lepton original one. 
In our model we have shown how the gravitational in¬ 
teraction allows a production of an asymmetry between 
HRN and HLA. This asymmetry is transferred to ordi¬ 
nary light neutrinos via deca y pr ocesses of the HRN and 
HLA (for further details see 133 . HH ). This happened at 
T ~ M ~ 10 12 GeV. Note that the necessary condition 
for our model is Td > M. Thus, neutrino asymmetry is 
transferred to lepton sector through these decays. The 
lepton asymmetry can be converted in the baryon asym¬ 
metry during sphaleron era, as shown in 1.1391]. 


B. Leptogenesis in ln(i?) 


M 

T d 

To^ 

up" 

“up" 

up" 

“up" 

up" 

“up" 

up" 


TABLE I: HRN mass and decoupling temperature in GeV for 
n = —0.01 in ln(R) gravities. 


M(GeV) 



FIG. 6: HRN mass in Ini? model in function of Td- 


We, at the beginning, impose n equal to —0.01 and Td = 
10 16 GeV. In this way LAF becomes: 


V = 


5.9 • 10 23 GeV 2 

M 2 


(V.52) 


In Table U we observe the variation of the mass of HRN 
in function of the decoupling temperature Td- In Figure 
[G]we show M vs Td- 

Besides we exhibit results in agreement with the pre¬ 
vious f{R) model for n = — 10~ 4 , observing that for 
Td = 10 15 GeV we get for the HRN a mass of the order 
of 10 12 GeV in agreement with the see-saw mechanism 
prediction (IV.501) as discussed before (LAF behaviour is 
in Figure 0. Straightforwardly both f(R) models intro¬ 
duced (f{R) = R + aR 2 and ln(R)) in order to explain 
leptogenesis give the expected LAF for the same value 
of the Temperature and HRN mass, i.e., Td ~ 10 15 GeV 
and M ss 10 12 GeV (how it is possible to see comparing 
results exposed in Figure [5] and 0. This is another im¬ 
portant hint that a little modification of the GR action is 
able to give the correct explanation of phenomenon such 
as leptogenensis. 


Hereafter, we show the leptogenesis in the frame of the 
ln(R) model of gravity, originally proposed here in this 
context. We proceed as it has been done before. We 
compute the same LAF given by (IV.481) evaluating R 
from (IIV.31I) . i.e., 


VI. RECONSTRUCTING THE POTENTIAL OF 
A PRIMORDIAL SCALAR FIELD FROM 
BARYO-LEPTOGENESIS 


n( 5 + 19n + 22n 2 + 8 n 3 )6T 6 
(2 + n) 2 M 2 Mp 


(V.51) 


The extra gravitational degrees of freedom generated 
from a f(R) theory of gravity can be represented from 
additional scalar fields. Here we aim to find the correct 

























potential form of a primordial scalar field that could be 
the main actor of leptogenesis and other phenomenon like 
Dark Energy or Inflation. 

For this scope, let us discuss the conformal transfor¬ 
mation applied to a generic f(R). Given a generic f(R) 
action it is possible, by means of conformal transforma¬ 
tion, to rewrite it in the Einstein frame, i.e., Ricci scalar 
plus a minimally coupled scalar field [9, I40l - l42 |. 

Hereafeter we consider the conformal transformation 
e 2x acting on the metric g^ v 0, i.e. 


Substituting this expression in (IVI.57I) . the potential be¬ 
comes 

V = 2 "" la(1 _7 rt) e fc ^[sinhM^ T -(VI.60) 
(an) "- 1 

In the case of kip -C 1 it is possible, by means of Tay¬ 
lor expansion, to assume a power law behaviour for the 
potential 0], i.e. 


9^ = e 2x <w . (VI.53) 


A particular choice of \ allow us to transform the begin¬ 
ning action f(R) in a new one composed by Ricci scalar 
plus minimally coupled scalar field [43l445l |. In particular 
the right choice of x i s: 


X = l -\n\f'(R)\. (VI. 54) 

Imposing 

kip = x, k = -^= , (VI.55) 


V ~ Vo<p 5 , (VI.61) 


with 


1 

2"- 1 a(l — n) r n 

Vo = — - —r^nr— , o =-- . 

(an)"- 1 n- 1 

In our case, i.e., n = 2 the potential (IVI.61I) assumes the 
following form: 



(VI.62) 


the Lagrangian density of f(R) can be rewritten in the 
(conformally) equivalent form [9|, RgI 

V=gW) = VH (- 5 # + - vj , 

(VI.56) 

where the potential V is defined as 

V = ■ (VI.57) 

Thus, we can explicitly compute the form of the po¬ 
tential V in the case of / = R + aR n . From (IVI.54|I . one 
gets 

/' = e 2 ^ , (VI.58) 


and knowing that f = 1 + anR n 1 , we can compute the 
Ricci in function of the new scalar field ip: 


R = 


— (e 2ktfi 
an 


i) 


(VI.59) 


where Vo = — 7 ^. While, always for n = 2 the general 
expression for the potential (IVI.60I) is 

— = e~ 2k<p sinh 2 (kip), (VI.63) 

Vo 

where Vo = — , 7= . The behaviour of such potential is 

4v2 Q; 

in Figure [5] 

A more interesting potential form is obtained for the 
f(R) = R + 7 R~ n (lnj^) m . In the case of m = 1 the 
Ricci in function of the scalar field is given by: 

/ nW ^». 1 /" (n+ i)(^-i)(^ e v)- j \ ^ 

(n + 1) (e 2kip — 1) 

V 

(VIM) 

where W(z) denotes the function that gives the principal 
solution for W in z = We w . The analytical expression of 


n—1 





























fi= 1 , y=l, m=l 


VII. ANALYSIS OF THE RESULTS AND 
CONCLUSIONS 
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V(v) 



FIG. 9: V vs ip for f(R ) = R+yR n (ln-S I ) rn where m, 7 , fi = 



FIG. 10: V vs for f(R) = R+'yR n {ln-S I ) m where m, 7 , ^ = 
1 , n = 10 -4 . 


the potential V{ip) becomes too complex (see Appendix 
A) and we report its behaviour for different values of 
the parameter n in Figure [9] We fixed the other pa¬ 
rameters, i.e., /x, 7 , because only n is responsible for the 
determination of the scale factor of the Cosmos as we 
have shown in Eq. (1IV.32I) . In Figure ITOl we represent 
the potential behaviour, for an order of magnitude of n 
equal to the one previously used, in order to explain the 
baryo-leptogenesis, i.e. n ~ 10~ 4 . Proceeding in this 
way, we can claim that the f(R) presents in Eq. (1IV.31I) . 
for to = 1 and n ~ 10 “ 4 , is conformally equivalent to 
a scalar-tensor theory of gravity where the scalar field, 
minimally coupled to the Ricci scalar, has a potential 
described by Figure [TUI Finally, it is worthwhile to point 
out that the general form of the two f(R), i.e., (IIV.24I) 
and (IIV.31D , or their conformally equivalent scalar-tensor 
theory with the potential described in this section, are 
able to generate baryo-leptogenesis, and they could also 
get rid of the problem of late and early time acceleration 
of the Cosmos. 


In this paper we have analysed the long-standing prob¬ 
lem of the origin of the asymmetry between matter and 
antimatter in the context of f{R) theories. According to 
Occam’s razor the solution here proposed, is the solution 
that needs fewer new conditions to solve the problems. 
In the context of f(R) theories of gravity, we realize the 
correct BAF/LAF only in presence of 3 conditions: ”CP 
violation, CPT dynamical violation, B/L violating reac¬ 
tion” . In particular the lepton asymmetry is generated by 
an interaction between chiral fermion and gravity (IV.3811 . 
This interaction splits energetic levels of neutrinos and 
antineutrinos in an expanding Universe under the force 
of a f(R) theory. Furthermore sphaleron converts this 
lepton asymmetry to the baryon sector. Besides, we have 
also described baryogenesis in an alternative way, in the 
context of an interaction between baryon and Ricci scalar 
[7s| which reproduces the expected BAF as consequence 
of a little modification of gravity a In this work we 
have proposed two different f(R), i.e. 

f{R) = R + aR 2 , (VII. 65) 

/ D \ m 

f(R) = R + 7 R~ n (in ^2 ) > (VII.66) 

where, in particular, the second one has never been stud¬ 
ied before in the context of the matter-antimatter asym¬ 
metry problem. 

Both model satisfy the conditions enumerated here¬ 
after, i.e. 

1 . e -C 1 in order to not change the standard thermal 
history of the Universe. 

2. reproduce LAF and BAF consistent with experi¬ 
mental data. 

Besides we have found the potential of a primordial 
scalar field, trough LAF constraints. Actually, f(R) 
gravities are conformally equivalent to a theory with tra¬ 
ditional Einstein term plus a scalar field. It is possible to 
find the potential, for the scalar-tensor theory equivalent 
to the f(R), trough conformal transformation of the met¬ 
ric. For both f(R) analysed in this paper we have found 
the potential of the scalar field generating the asymmetry 
between matter and antimatter. 

It is worthwhile to highlight how f(R) theories of grav¬ 
ity can introduce a small perturbation to the GR scale 
factor that allows us to obtain the expected value for the 
asymmetry factor. 

At the end we point out that, as shown in this paper, 
f(R) theories may be the ultimate solution for most of 
open problems in modern cosmology, e.g. Dark Energy, 
Inflation, Dark Matter, Bario-Leptogenesis. 
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Appendix A 

The general expression for the potential of the scalar field action conformally equivalent to f(R) = R + 'yR~ n ln(-^) 


is: 


where 


I'M = 


(.67) 


B = 2 


f 

^jnW ^ 

'(„+l)(e 2 ^- 1 )(^e 1 /") B+1 \ \ 

n+l \ 


7 n y 




(n + 1) (e 2kv — 1) 


V 

V 

) 

) 


and 


A = 


(n+ 1) log 


— qn log 

V 

\ 


■ynW ^ 

(n+ 1 )(e 2 '=*’- 1 )(^2 e l/»)’ l + 1 \ \ 

n + l 

7 " J 1 


(n+ l)(e 2fe( ^-l) J 


d 2 


V V 


(n + l)(e 2/e <^ -l) 


) 


\ 


- 1 


) ) 


( jnW ( ^ 2kv -^ 1/n r +1 


(n+ l)(e 2fe v=-l) 


7 nW 


(V 


(n+l)(e 2fc ^-l) 


\ n+lN l 

\ 

) 1 

/ 


(. 68 ) 


(.69) 


W(z) is the function that gives the principal solution for IT in z = and in Wolfram Mathematica® is named 

ProductLog[z\. _ 



